Contributed paper to Lepton Photon 2001 
Rome, July 23-28. 

Detection of ^ /o(980)7, (/> -^ ao(980)7 into 5 photons with 

_ KLOE at DA<I>NE 

o 
o 

C^ . The KLOE Collaboration 



o 



> 



X 



^ 



A. Aloisio? F. Ambrosino^ A. Antonellif M. Antonellif C. Bacci," G. Barbiellini? F. Bellini," 
G. Bencivenni,'^ S. Bertolucci,'^ C. Bini,'' C. Bloise,"^ V. Bocci,''' F. Bossi,"^ P. Branchini," S. A. Bulychjov,^ 
G. Cabibbo,*-' R. Caloi,'^ P. Campana,'' G. Capon,'^ G. Carboni™ M. Casarsa,' V. Casavola,^ G. Cataldi,^ 

F. Ceradini," F. Cervelli^ F. Cevenim,^ G. Chiefari,^ P. Ciambrone,'^ S. Conetti,'' E. De Lucia,*^ 

G. De Robertis," P. De Simone,^ G. De Zorzi,*-' S. Dell'Agnello,'= A. Denig,'= A. Di Domenico,'' C. Di Donato,^ 
(;S1 ■ S. Di Falco,'' A. Doria,^ M. Dreucci,'= O. Erriquez," A. Farilla," G. Felici,^ A. Ferrari," M. L. Ferrer,^ 
^D I G. Finocchiaro,'^ C. Forti,"^ A. Franceschi,'^ P. Franzini ,*"''* C. Gatti/ P. Gauzzi,*^ A. Giannasi/ S. Giovannella,'^ 
P: ■ E. Gorini,'^ F. Grancagnolo,'^ E. Graziani," S. W. Han,^^'' M. Incagli^ L. Ingrosso,^ W. Kluge,'* C. Kuo,'' 

V. Kulikov,-'' F. Lacava,*^ G. Lanfranchi,'^ J. Lee-Franzini,'^'° D. Leone,'^ F. Lu,'^"'' M. Martemianov,'^'-'' 
M. Matsyuk,^'/ W. Mei,^ A. Menicucci';" L. Merola,^ R. MessiT* S. Miscetti,^ M. Moulson,'^ S. Miiller,'' 
F. Murtas,"^ M. Napolitano,^ A. Nedosekin,'^'^ M. Palutan," L. PaoluziT* E. Pasqualucci,''' L. Passalacqua,'^ 
D I A. Passeri," V. Patera,'''^ E. Petrolo,'' D. Picca,'' G. Pirozzi,^ L. Pontecorvo,'^ M. Primavera,*^ F. Ruggieri," 

O ' P. Santangelo,'^ E. Santovetti™ G. Saracino,^ R. D. Schamberger," B. Sciascia,*^ A. Sciubba,''^ F. Scuri,'' 

<D '. I. Sfiligoi,'^ J. Shan,^ P. Silano,*^ T. Spadaro,'= E. Spiriti," G. L. Tong,^'^ L. Tortora," E. Valente,'^ P. Valente,'= 

^ B. Valeriani,'' G. Venanzoni^ S. Veneziano,*^ A. Ventura,'^ Y. Wu,'=^'' G. Xu,'^''' G. W. Yu,'=^'' P. F. Zema/ 

Y. Zhou'^ 

/\ ' " Dipartimento di Fisica deU'Universita e Sezione INFN, Bari, Italy. 

^ , ^ Permanent address: Institute of High Energy Physics of Academica Sinica, Beijing, China. 

'^ Laboratori Nazionali di Frascati dell'INFN, Frascati, Italy. 
'^ Institut fiir Experimentelle Kernphysik, Universitat Karlsruhe, Germany. 
'^ Dipartimento di Fisica deU'Universita e Sezione INFN, Lecce, Italy. 

f Permanent address: Institute for Theoretical and Experimental Physics, Moscow, Russia. 
^ Dipartimento di Scienze Fisiche deU'Universita e Sezione INFN, Napoli, Italy. 
' Physics Department, Columbia University, New York, USA. 
^ Dipartimento di Fisica deU'Universita e Sezione INFN, Pisa, Italy. 
'^ Dipartimento di Fisica deU'Universita "La Sapienza" e Sezione INFN, Roma, Italy 
' Dipartimento di Energetica deU'Universita "La Sapienza" , Roma, Italy. 
™ Dipartimento di Fisica deU'Universita "Tor Vergata" e Sezione INFN, Roma, Italy 
" Dipartimento di Fisica deU'Universita "Roma Tre" e Sezione INFN, Roma, Italy 
° Physics Department, State University of New York at Stony Brook, USA. 
■^ Dipartimento di Fisica deU'Universita e Sezione INFN, Trieste, Italy. 
^ Physics Department, University of Virginia, USA. 

Abstract 

The KLOE detector at DA$NE, the Frascati (^-factory, has collected about 30 pb~^ by the end of 
2000, corresponding to about 90 millions of <^(1020) mesons produced. The five photon final state has 
been exploited to study the rare decays of (j> into /o (980)7 s-nd ao (980)7, with statistical accuracy never 
reached before. The preliminary results for these branching ratios, from a data sample of 17 pb~^. 



are the following: Br((j!> -> /o(980)7) = (23.7±0.6«tat)x 10-^ Br{4> -^ ao(980)7) = (5.8±0.5^tat)xl0-^ 
The systematic error is expected to be less than 10%, a precise evaluation is going on. For the ratio 
of the branching ratios of (;A — > /07 to ^ ao7 we find 4.1±0.4stat- 

1 Introduction 

In the last years there has been a revival of interest in the lowest-lying scalar mesons with masses below 

1 GcV which has motivated the experimental effort going on with the KLOE detector |]l| at DA$NE|p| 
the Frascati (/)-factory, where these resonances can be studied in detail in (f> radiative decays. 

The 1 GeV region is an interesting and challenging domain: on one side it is far below the perturba- 
tive QCD regime, on the other side strict Chiral Perturbation Theory is not expected to make reliable 
predictions at these energy values where resonance effects can show up. This region has thus become 
an important test ground both for analyses based both on QCD sum-rules[B| and for effective chiral 
modelsgl, I. 

The controversial nature of the /o(980), ao (980)0 and the poor knowledge of their properties adds 
further interest to this 1 GeV energy region: along the years several proposal have been suggested 
concerning the constitution of these scalars as complex qqqq states [g|, KK molecules § or ordinary qq 
mesons [|lO| . A precise measurement of Br(0 -^ fol), Br((/) -^ 097) and of the ratio of these branching 
ratios can help in better establishing the nature of these mesons [|ll[. 

In this paper we describe the preliminary results on the decays </> — > /07, 097 — ^67 from the analysis 
of 17 pb~^, corresponding to about 50 millions of </> decays, from the data sample collected by KLOE in 
the year 2000. A peculiar aspect of these decays is that both the branching ratio and the resonance mass 
shape depend strongly on the meson structure. Therefore, the analysis follows a scheme independent as 
much as possible on the model implemented in the Montecarlo. 

2 Selection criteria for five photon events 

The events are characterized by the presence of five prompt photons, i.e. photons coming from the 
interaction point (LP.) of DA<I>NE. These photons are detected as energy deposits in the calorimeter that 
satisfy the condition t — r/c = 0, where t is the arrival time, r is the distance of the cluster from the LP. 
and c is the speed of light. We define a photon to be "prompt" if \t — r/c\ < 5crt[|l2|. An acceptance 
angular region corresponding to the polar angle interval 21°— 159° is defined for the prompt photons, in 
order to exclude the blind region around the beam-pipe. The processes that mainly contribute to the five 
photon final state are the following: 

- </> -> /o7 -^ 7r'^7r°7 

- ^ p"^o ^ 7r"7r07 

- (f> -^ ao7 -^ riTT°-f (?7 -^ 77) 

- (j> ^ P^t:^ — > ?77r*'7 (77 -^ 77) 

and the non resonant ones: 

- e+e~ -^ LOTT^ — ^ 7r"7r"7 

- e"'"e~ — > lot:'^ — > 777r''7 {-q — > 77) 

Also the three and seven photon final states 

- (j) —>■ rjj ^ 777 



- — !■ TT^'y — > 777 

- e+e" -^ 77(7) 

- (j> ^ rjj —^ vr^vr*' 71*^7 

have to be taken into account, because they can simulate five photon events due to machine background 
and cluster splitting or merging of close photons and loss of soft photons. 

The analyses described in this paper make use of a constrained fit ensuring kinematic closure of the 
events. The free parameters of the fit are: the three coordinates (x, y, z) of the impact point on the 
calorimeter, the energy and the time of flight for each photon coming from the LP., the two energies of 
the beams and the three coordinates of the position of the LP.. The analysis procedure adopted is the 
following: 

1. fully neutral events with exactly five prompt photons are selected; 

2. the kinematic fit is applied on these events a first time with the constraints of the total energy and 
momentum conservation and satisfying t — r/c — for each prompt photon; 

3. the photon pairing is performed looking for the best combination, by minimizing the appropriate 
X^ for each of the following hypotheses: tt^tt'^j, V'^'^li ujtt'^, rij — > 777 and tt'^j. 

4. the kinematic fit is applied again to each of the best combination constraining also the invariant 
masses of the photon pairs assigned to tt*^ and ry, without any assumption on the /o and oq mass. 

3 ^ 7r^7r^7 

The main background channels, with the expected signal to background (S/B) ratio obtained from |^, 
|l4| , [l5| , [lq| , are listed in Tab. ^ After applying the first kinematic fit, the best photon combination 
producing the two 77 pairs with 7r° mass is searched for with the additional condition M-^,^ > 700 McV 
(the expected /o mass region). According to our simulation, this procedure correctly assigns photons in 
95% of the cases. The analysis cuts to select 4> -^ /o7 ^ 7r"7r°7 events are: 

1. xL^/ndf < 5; 

2. 5(7 cut on the reconstructed tt^'s mass (the output of the first kinematic fit, without mass con- 
straints, is used); 

3. cosV'ojir > 0.4, where ^ is the angle between the primary photon and the pion's flight direction in 
the tt'^tt^ rest frame. This cut is performed to reduce the —^ ao7 and — > p^n^ backgrounds; 

4. x^Tr/'^df < 3 and a 3cr cut on the reconstructed u mass to veto e^e^ — > wtt'^; 

5. Etot > 900 MeV and 175° < A0 < 185° for the two most energetic photons to reject e'^e^ -^ 77(7); 

6. 3 cr cut on the radiative photon energy {E^ rad ~ 500 MeV) to reject (f> — > tt'^j. 

In order to evaluate the analysis efficiency as a function of M^^, the generated /o mass has been 
divided into 20 MeV bins. The total variation of the efficiency does not exceed 30% of the mean value, 
which is reported in tab. M together with the efficiencies for the background channels. 

A fit procedure is in progress to evaluate the branching ratio and to obtain the /o mass shape by 
folding the available theoretical dN/dM shapes with the reconstruction effects. Since this procedure is 
not completed and a reasonable /o mass shape was chosen in the Montecarlo, the branching ratio is 
evaluated using bin by bin efficiency while data and Montecarlo distributions are compared as follows: 
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Figure 1: — > /07 — > 7r'^7r^7 distributions: x^/ndf (up-left), tt'^tt'^ invariant mass (up-right), polar 
angle (down-left) and cost/i (down-right). Black circles are data, solid histograms are the expected 
MC signal-t-background spectra while the background contributions are displayed in colors: green is 
e'^e^ — » ciJTr*^, yellow is cj) —>■ rjj, magenta is cj) —^ pn^ and blue is — > 007. The MC /o shapes have been 
obtained by weighting the Mj^^^ bin contents to reproduce the shape of the data. 
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Figure 2: (j> ^ /07 -^ tt^tt'^j events: energy distributions (up) [left: radiative photon, right: 7's from 
pions] and the two 7r° masses (down). Black circles are data after background subtraction, solid histograms 
are the expected spectra for the signal. The MC /o shapes have been obtained by weighting the M-^^^ bin 
contents to reproduce the shape of the data. 



Decay channel 


Natural S/B 


Analysis efficiency 


^ /o7 ^ 7r°7r°7 


— 


39.7% 


e+e- ^ Lun" ^ 7r07r07 


0.6 


1.2% 


^ p"7rO ^ ^"7r07 


3.7 


4.9% 


^ -^ ao7 ^ 777r*'7 ^ 777r'^7 


3.5 


1.9% 


(j) ^> rj"f ^ 7r"7r''7r''7 


0.02 


5 X 10--* 



Table 1: S/B ratios and analysis efficiencies for the (f) -^ foj -^ tt^tt'^j decay and related background. 



1. 



/o7 MC events are analyzed in slices of Mt^tt to obtain the various histograms shapes; 



2. after subtracting the expected background from the observed Af^^ distribution, for each M^^ bin 
we evaluate the scale factor that, apphed on the MC M^^r shape, reproduces data; 

3. these scale factors are then applied to all MC "sliced" histograms which are then summed to produce 
the MC distributions to be compared with the data. 

This implies that any effect of bad event reconstruction that gives rise to a distorted M^^ shape is 
neglected. We expect such contribution mainly from bad photons' reconstruction ('^ 4 -^ 5%) and wrong 
7's pairing. The data-Montecarlo comparison for the most relevant distributions after this procedure 
is shown in Fig. |l|. The behaviour of the energy spectra and of the pion masses after the background 
subtraction is also reported (Fig. g). 

Applying the analysis to the 17 pb^^, 1967 events survive while the Montecarlo expected background 
contribution is 305± 13stat ■ The corresponding branching ratio has been evaluated by normalization with 
respect to the (/> — » 777 ^ 777 events, using the PDG |16| value for the Br((/) — > 777 -^ 777)- Neglecting 
the interference between the signal and (f> —i- pir'^ — > tt'^tt'^j, the final branching ratio is: 



Br( 



/o7 ^ ^"^"7) = ( 7.9 ± 0.2,tat ) X lo- 



ci) 



for AU^ > 700 MeV. 

The systematic error is under study and it should not exceed 10%. 



4 ^ ?77r^7 

The (f) — > 7771*^7 events, with 7/ — > 77, have been considered to study the cf) — > ao7 decay. The corresponding 
S/B ratio of the signal and the main background channels are reported in tab|2[ 

A first selection to reject e"''e~ — > 77(7) events has been applied. Then the 777 — > 777 contamination 
is eliminated by rejecting events with a two photon invariant mass near the rj mass together with another 
photon energy near 363 MeV. 

In order to reduce the n'^TT^'y background, the photon combination obtained for the e"'"e~ — > ujtt'^ — > 
7r°7r°7 hypothesis has been exploited. The correlation between the angle tp (angle between the non asso- 
ciated photon and the primary tt^ in the dipion rest frame), and the energy of the non associated photon 
is shown in fig.0.a,b for MC samples and in fig.H.c for data. From the comparison several contributions 
can be identified. In order to avoid introducing any bias in the unknown oq shape, the analysis cuts 
have been studied using a (p —> ijir^'j MC sample with an almost flat 777r*' invariant mass. That sample is 
compared in fig.0.d with the dominant background; also the applied cut is shown. 

Since the cj) — > /5°7r° — > 7r°7r°7 events populate the same region of the scatter plot of figJq as the signal, 
this contamination should be reduced by using different variables. The choosen ones are shown in the 
scatter plot of fig.0: the difference between the x^ of the second fit in the hypotheses 7r'^7r°7 and 7777^7 
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Figure 3: V' vs E-y - (a) Data; (b) MC: black - ojtt'^ — > tt'^tt'^j, red - ao7 -^ w'^li green - fo7 -^ 7r°7r''7, 
blue - rj'j ^ 777; (c) MC: black - W7r° -^ 7r°7r°7, red - p°7r° -^ 7r'^7r"7, green - p"7r° -^ iT^^l^ blue - 
uiTT^ -^ ?77r'^7; (d) MC: black - ujtt^ -^ tt'^tt'^j; red - rjir'^'-f flat (see text); the accepted region is below the 
solid line. 
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Figure 4: t] reconstructed mass vs xi-rr^y ~ X^-n-t for MC samples: red - -q-K^^ flat; blue - p°7r° -^ 7r"7r°7; 
the applied cut is also shown. 
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Figure 5: (a) Reconstructed 77 mass divided by the resolution; full histogram: 777 -^ 7r°7r'^7r'^7 (MC); (b) 
Reconstructed tt*^ mass divided by the resolution. 
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Figure 6: (a) jyTr" invariant mass - points: data, histograms: background spectra (MC), red - tt'^tt^j, 
green - e"'"e~ — > wtt" — > rjn'^j, blue - rij —> tt'^tt'^it'^j, black - total; (b) rjir'^ invariant mass background 
subtracted; (c) cos^ distribution - points: data, histograms: background, signal plus background. 



Process 


Natural S/B 


Final efRciency 


-^ ao7 -^ V^^l 
e+e- — * W7r° -^ ijtt^j 


5.3 

71 


27.2% 
27.1% 
25.8% 


e+e~ — > a;7r° — + tt^tt^j 
-^ /o7 ^ 7r°7r°7 


0.14 

1 
0.27 


3.5x10-3 

5.0% 

1.4-3.0xl0~3(*) 


-^ 777 ^ 777 


6.1x10-3 
7.5x10-3 


3.3xl0-«^ 
5.4x10-* 



Table 2: S/B ratios and efficiencies for i-jit'^j decay and related background. 
*^*)Depending on the /o shape. 



and the reconstructed rj mass. The accepted region is defined by x ~ X >3 plus the internal 
part of the ellipse. 

Last contamination that has to be reduced is — > 777 — > tt^tt^tt'^'^. In fig.@.a is reported the invariant 
mass of the two photons associated to the ij; the comparison with the MC shows that the background 
events mainly populate the tails of the distribution. A two standard deviation cut is then applied to 
reduce the 777 contamination. The final efficiencies are reported in tabM 

At the end of the selection procedure 666 events survive. The expected number of background events 



is 253±11 mainly due to 



pO^o 



7r*'7r*'7 and 



VI 



^OttOtt^. 



The invariant mass of 7777*^ is 



shown in fig.o, together with the cost? distribution for the non associated photon. The corresponding 
branching ratio is: 

Bricj) ^ 777r%) = (7.4 ± 0.5stat) x 10-^ (2) 

The systematic error is still under study; according to preliminary evaluations it should not exceed 10%. 
This result is in agreement, within the errors, with the values published by the VEPP-2M experiments ]l^, 
^ and with the KLOE resuh of the analysis of the 1999 statistics^. 

The branching ratio (g) includes the cj) — » p'^ir^ -^ 7771*^7 contribution. If the interference of this process 
with 4> ^^ ao7 ~^ 7771*^7 is negligible [|l9| , one can subtract the p'^n'^ contribution. To evaluate the number 
of expected p'^n'^ events, the average of the very recent measurements of the Br(73° -^ 777) performed 
by the two VEPP-2M experiments pof has been used. This average value turns out to be Br(p° -^ 
777)=(3.0±0.3)xl0-* (which is different from the PDG value(l|l Br(p° -> 777)=(2.4±0.9)xl0-*), from 
which follows an expected number of 86±9 p^n'^ -^ rjir^j events. After the subtraction one obtains: 



Br{(f> -^ ao7 ^ V^"l) = (5-8 ± 0.5,tat) x 10-^ 



(3) 



/o7 to 



ao7 ratio 



From the results described in the previous sections a preliminary estimate of the ratio of the two branching 
ratios can be obtained. Assuming that 7r*'7r*'7 and 7771*^7 are the dominant decay modes of respectively /o 
and flQ and using the values of eqs.(P and (|^), one obtains (with Br((/) — > /o7)=3Br(0 — > /07 -^ n^n^^)): 



Br{4> -> /07) 



Br{(l) -^ ao7) 
in reasonable agreement with the prediction given in Jll| 



= 4.1 ±0.4, 



(4) 



10 



6 Conclusions 

Events with five prompt photon have been selected from a sample of 17 pb^^ of the data collected by 
KLOE in 2000, corresponding to about 50 millions of (p meson produced. From the analysis of the Tr'^n'^'y 
channel, a Br((/) -^ /o(980)7 -^ tt'^7t'^^)~{7 .9±0.2stat)'x 10~^ has been obtained. The result of the analysis 
of the r/Tr'-'j with 77 ^ 77 channel, is: Bi-{(p -^ ?77r°7)=(7.4±0.5stat)^10"^- By subtracting the expected 
(/) -^ p°n° — > ?77r°7 contribution, one can evaluate Bt{(/) -^ ao(980)7 -^ r]n°^)— {5. 8±0. 5 stat)x^O^^- The 
ratio of the branching ratio (/) — > /07 to — > ao7 = A.liOAstat is obtained. 
Work is in progress to fit the spectra and to extract the /o and oq resonance parameters. 
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